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Abstract. In a previous work, we used a polarization condition to show 
that there is a family of U(l) symmetric solutions of the vacuum Einstein 
equations on E x S 1 x R (E any 2 dimensional manifold) such that each 
exhibits AVTD 1 behavior in the neighborhood of its singularity. Here we 
consider the general case of S 1 bundles over the base Exi? and determine a 
condition, called the half polarization condition, necessary and sufficient in 
our context, for AVTD behavior near the singularity. 

1 Introduction. 

A rigorous study of the singularities in cosmological solutions of the vacuum 
Einstein equations has been hampered by the fact that the generic such 
solution is expected to have a singularity of the oscillatory type predicted by 
Belinsky, Lifshitz and Khalatnikov [BLK]. There is currently no satisfactory 
mathematical method for treating such oscillating singularities, at least when 
the spacetimes under study are spatially inhomogeneous. For this reason 
much effort has gone into the study of families of solutions which have milder 
cosmological singularities such as those of AVTD (asymptotically velocity 
term dominated) type, for which rigorous, so called Fuchsian, methods are 
available. 

To suppress the oscillatory behavior expected for the generic solution, 
one can 
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• (i) introduce suitable matter sources such as scalar fields and study the 
solutions of the associated non vacuum field equations [AR], 

• (ii) study higher dimensional models motivated by string or super- 
gravity theories wherein (for sufficiently high dimensions at least) the 
oscillations are naturally suppressed [DHRW]; or 

• (iii) remain in 3+1 dimensions but impose a combination of symme- 
try and polarization conditions in order to achieve the desired AVTD 
behavior. 

For the case of U(l) symmetric vacuum solutions on the trivial S 1 bun- 
dle T 2 x R x S* 1 — > T 2 x R (with U(l) symmetry imposed on the circular 
fibers) Isenberg and Moncrief [IM] have showed, using Fuchsian methods, 
that AVTD behavior is achieved provided the solutions considered are at 
least half polarized in a certain well defined sense. The half polarization con- 
dition includes, as a special case, the fully polarized solutions wherein the 3 
planes orthogonal to orbits of the U(l) isometry action are integrable and the 
vacuum 3+1 field equations reduce to a system of 2+1 dimensional Einstein 
equations coupled to a massless scalar field on the quotient manifold T 2 x R. 
The more general (half polarized) solutions admit, in addition, half the extra 
(asymptotic) Cauchy data expected for a fully general, non polarized solu- 
tion of the same (U(l) symmetric) type. On the basis of numerical studies 
due to Berger and Moncrief the fully general, non polarized U(l) symmetric 
vacuum solution on T 2 x R x S* 1 — > T 2 x R is expected to have an oscillatory 
singularity and hence not to be amenable to Fuchsian analysis [BM]. 

Choquet-Bruhat, Isenberg and Moncrief have extended the analysis given 
in [IM] to cover the case of polarized U(l) symmetric vacuum solutions on 
manifolds of the more general type S 2 x R x S 1 — > X 2 x R, where E 2 is an 
arbitrary compact surface and the bundle (in view of the assumed polariza- 
tion condition) is necessarily trivial. In the present paper the polarization 
restriction is eliminated in favor of an appropriate half polarization condition 
and the limitation to trivial S* 1 bundles over the base S 2 x R is also removed. 
The present work thus demontrates the existence of a large family of vac- 
uum U(l) symmetric solutions of half polarized type defined on trivial and 
non trivial bundles over X 2 x R (with S 2 an arbitrary compact surface) and 
having AVTD singularity behavior. The half polarization condition used in 
[IM] involved requiring one of the asymptotic functions to vanish. The half 
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polarization condition which we find here necessary and sufficient for possi- 
ble AVTD behavior can be understood in terms of the behavior of the VTD 
solutions to which our solutions converge as one approaches the singularity. 
Specifically a VTD solution is half polarized if and only if the set of geodesies 
in the Poincare plane which represent it (at different spatial points) all tend 
to the same point as t appoaches the singularity. 

2 Einstein equations. 

A spacetime metric on a manifold V4 = M x R, with M an S* 1 principal fiber 
bundle over a surface E, reads, if it is invariant under the S 1 action on V4, 

^g = e~^ Wg + e^idO + A) 2 (2.1) 

with 6 a parameter on the (spacelike) circular orbit, <fi a scalar, A a locally 
defined 1 - form and ^g a lorentzian metric, all on V3 := E x R. 

The vacuum 3+1 Einstein equations Ricci(^g) = for such an S* 1 sym- 
metric metric on V 4 are known [M86], [CB-M96] to be equivalent 2 to the 
wave map equation from (V3/ 3 ) g) into the Poincare plane P =: (R 2 ,G), 
$ = ((f), uj) : V 3 — > R 2 , where 

G = 2(d<f)) 2 + -e-^(duj) 2 , (2.2) 

coupled to the 2+1 Einstein equations for ^g on V3 with the wave map as 
the source field. The scalar function u on V 3 is linked to the differential F 
of A by the relation 

duj = e 4<t> *F, with F = dA. (2.3) 

Thus in local coordinates x a , a — 0, 1, 2, on V3, with rj the volume form of 

(3)g _(3) g a pdx a dxP , one has 

F aP = ^e-^ x d X io. (2.4) 

The wave map equations are, with ( 3 ^V the covariant derivative in the 
metric ^g 

g^i^VM + l -e-^d a udpuj) = (2.5) 

2 If we choose an arbitrary harmonic 1 - form appearing in the solution to be zero. 
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g a/3 { {3) V a d^ - Ad a ud P 4>) = 0. (2.6) 

The 2 + 1 Einstein equations are, with "." indicating a scalar product in 
the metric G 

^R af3 = d a <Z>.dp<Z>. (2.7) 

To solve these equations we choose for a zero shift, we denote the 
lapse by e A and we weigh by e A , without restricting the generality, the t 
dependent space metric g = g ab dx a dx b , a, b = 1, 2. That is, we set 

^g=-N 2 dt 2 + g ab dx a dx b with N = e\ g ab = e x a ab , (2.8) 

We denote by a ab the contravariant form of a. The extrinsic curvature of E t 
in (V 3 &g) is 

1 1 

Kb ■= -yjydtgab = -^abd t \ + d t a ab ). (2.9) 
The mean extrinsic curvature r is therefore 

r:=g ab k ab = -e~\d t \ + ^), (2.10) 

where we have defined 

V := a ab d t a ab . (2.11) 

The connection coefficients (Christoffel symbols) of ^g are found to be 
(note that ( 3 V° = -e~ 2A , ^g ab = g ab = e' x a ab ) 

{3) K b = Kb(g) = r°M + ^( W + W - a^^A) (2.12) 
wrgo = ftA, (3) r° a = 9 a A, ( 3 )r- = ^e A 9 a A, (2.13) 

(3) r°, = ~e- x k ab , ^T b a0 = -eX- (2-14) 
In particular it holds that 

( V W) C = \\>e~ 2 \ (2.15) 
We see that the metric ^g is in harmonic time gauge if and only if ip — 0. 
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The Einstein equations split into constraints and evolution equations. We 
denote by S% =^ Rg — \5p^R the Einstein tensor of ^g, by Tp the stress 
energy tensor of <J>, and we set E^j = Sp — Tp. The constraints are: 

C = TP = - l -{R{g) - k.k + r 2 - e~ 2X d t <$>.d t <$> ~ g ab d a $.d b $} = (2.16) 

and (indices raised with g ab , V the covariant derivative in the metric g) 

C a = e A E° = -{V b k b a - d a r + e - A <9 a $.<9 t $} = 0. (2.17) 

The evolution equations are, with iV = e A , 

N{ i3) R b a - p b a ) ee -d t k b a + Nrk b a - V b d a N + NR b a - Nd a <5>.d b <5> = 0. (2.18) 

In order to obtain a first order system in the Fuchsian analysis that we 
will make, we introduce auxiliary unknowns $t, $ a , a ab which are identified 
with the first partial derivatives of $ and the covariant derivative of a with 
respect to a given t independent metric a. These new unknowns satisfy the 
evolution equations 

dt$ = $ t , (2.19) 

d t $ a = dj&t, (2.20) 

dtof = V c d t v ab (2.21) 
where, by the definitions of a and k, 

d t a ab = 2e 2X k ab + a ab d t X. (2.22) 

The function A is not left unknown, but rather is determined by a gauge 
condition from its VTD value. 

3 VTD equations and solutions. 

The Velocity Terms Dominated equations are obtained by dropping the space 
derivatives in the equations. 

We denote by a tilde quantities which are independent of t, and we denote 
VTD solutions using a hat. 
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3.1 Einstein evolution VTD solutions. 

In order to obtain a global (on E) formulation we choose a VTD metric which 
remains in a fixed conformal class over E as t evolves; we set 

bob = Oab and g ab = e x a ab . (3.1) 

Then 

V> = 0, g ab = e" V b , d t g ab = e X a ab d t \, (3.2) 
and the definition of k gives that 

k ab = -h ab dX k b a = -h-H%\ f:=k a a = -e-' x d t X. (3.3) 

Requiring that these VTD quantities satisfy the VTD evolution equations, 
we obtain 

d t k b a = Nrk b a . (3.4) 
Therefore, by straightforward computation 

d 2 t X = 0; hence A = A — vt (3.5) 

with A and v arbitrary functions on E, independent of t. Then we have 

k ab = \va ab , k b a = hje-Hl r = e-h. (3.6) 

3.2 Wave map VTD solutions. 

The results for a VTD wave map are very different from the results obtained 
for a scalar function [CBIM]. If we drop space derivatives in the wave map 
equations we obtain geodesic equations in the target manifold, with t the 
length parameter on these geodesies so long as the 2+1 metric is in harmonic 
time gauge. If we make the change of coordinates Y = e 2 ^ in the target 
(which defines a diffeomorphism from R 2 onto the upper half plane Y > 0), 
the metric G takes a standard form for the metric of a Poincare half plane; 
namely 

G ,i { ^! } , y = e * (3 . 7) 

The VTD, geodesic, equations written in this metric read, with a prime 
denoting the derivative with respect to t 

u" - 2Y- 1 u'Y' = 0, (3.8) 
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Y" + Y^uo'X' = 0. (3.9) 

The general solution of these geodesic equations is represented in these co- 
ordinates, as is well known, by half circles 3 centered on the line Y = 0; 
specifically, with A and B arbitrary constants (that is, independent of t) , the 
solution takes the form 

Co = B + Acos6, Y = Asm6, O<0<tt. (3.10) 

These functions uo and Y satisfy the differential equations 3.8, 3.9 if and only 
if it holds that: 

¥-^e e - (3 - n) 

Integrating this equation we have that, with w independent of t 

d' = -wsm6. (3.12) 
Another integration gives that, with independent of t 

tan- = ee _iW . (3.13) 

If we now make the substitution A = e 2<t> and B = Co, then 3.10 reads 

= 0+ -log(sin9), Co = Co + e 2 ^ cosO . (3.14) 

Remark 3.1 The set of above formulas is identical to the following one 

Co = Co + e 2 ^ = , e 2(p = Y = e 2 ^ = . 3.15 

l + 2 e- 2 ^ l + 2 e- 2?5t 

Y tends to zero when t tends to oo, but Co tends to Co + e 2 ^. 



3 We discard here the special case which corresponds to the polarized case, treated 
elsewhere, where these circles are centered at infinity. The geodesies are then the half 
lines X= u =constant,. 
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3.3 Einstein constraint VTD solutions. 

We deduce from 3.14 and 3.12 that 

^0 = \y- 1 Y' = = --wcosO (3.16) 

Y 2 2 sin6 2 v ; 

and 

e - 2 h t Cj = -& = wsinO. (3.17) 
The Einstein VTD constraints reduce to the following equation 

2C = -k.k + r 2 - e- 2 ~ x {2(d t 4>) 2 + \e~ A \d t u) 2 } = 0. (3.18) 

We have, using 3.16 and 3.17, 

2(^0) 2 + l -e-^{d t u) 2 = l -w 2 . (3.19) 

We deduce therefore from 3.6 that the VTD constraint 3.18 is satisfied if and 
only if 

v 2 = w 2 (3.20) 

4 Fuchsian expansion. 
4.1 2+1 metric expansions. 

For the unknowns a and k we choose the following expansions, with the 
various e's being positive numbers to be chosen later 

a ab = ~ab + e -e a t u ab ^ 

kl = e-\f5 b a + e-^u{ a ). (4.2) 

Then 

r^k a a = e~\v + e-^ul a ). (4.3) 
We take as a gauge condition 

A = A; hence d t \ = —v. (4.4) 

Comparing the expressions 4.4 and 2.10 for r, we find that this condition is 
equivalent to the gauge fixing requirement 

e"%+^ = 0. (4.5) 
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4.2 Wave map expansion. 

We expand <£> near its VTD value; that is we set 

= + e~ e ^-u^ with = + ]-log(sin0), (4.6) 
while for u, for convenience of computation, we choose to set 

oj = Qj + e^e -6 "*^ with Cj = uj + e 2 ^cos6, (4.7) 

4.3 Expansion for first derivatives. 

We expand the auxiliary unknowns near the values of the derivatives of the 
VTD solution. That is we set (see 3.16, 3.17) 

4> t = d t 4> + e-^u^ = - l -wcosQ + e~ £ ^% t , (4.8) 

out = d t Cj + e 24> e- £ ^u wt = e 2 *wsin 2 6 + e^e'^u^ (4.9) 
The expansions of a and u a are defined similarly by setting 

0a = <9 a + e- £ *' U4>a , oo a = d a u + e^e^u^. (4.10) 

We next compute <9 a and d a u. It follows from 3.14 that 

d a <t> = d a (t>+ 1 —d a 6, d a u} = d a u + e 2 < t> (2cos6d a <P-sin6d a e). (4.11) 

We compute d a 9 using 3.13 and elementary properties of sine and cosine. We 
find that 

d a = Q- l sin0d a {Q - tw) (4.12) 

Therefore it holds that 

-s ~ ~ CO **>{/ ~ 

d a <p = d a <p + e- 1 — 9 (e - wt). (4.13) 

Then, writing d a Co as sum of a term independent of t plus terms tending to 
zero when t tends to infinity, we have 

d a tb = d a (Co + e 2 *) - e 2 ^[2(l - cos9)d a ^ - e _1 sm 2 ^ a (e - wt)] (4.14) 

For a^ b , since V c a ab = 0, we set 

of = e-^'*< c . (4.15) 
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5 Fuchsian system for the evolution equations. 



Given the Fuchsian expansions of the previous section, the Einstein - wave 
map evolution system reads as a first order system for the set of unknowns 
U = (u a , u k , u$ t , u a >). 

The differential system for U is Fuchsian in a neighbourhood of t — +00 
if it takes the form 

d t U-LU = e-' a F(t,x,U,dU) (5.1) 

with L a linear operator independent of t with non negative eigenvalues, \x 
a positive number and F a set of tensor fields linear in dll, continuous in 
t, analytic in x and U and uniformly Lipshitzian in all its arguments in a 
neighbourhood of U — 0, for t large enough. 

5.1 Einstein evolution equations. 

5.1.1 Equation for u a . 

The Fuchsian expansion 4.2 for k yields the following equation: 

d t g ab = 2Nk ab = 2e x g ac k b c = e-\vcx ab + 2 e - £fc V ac w 6 fcc ). (5.2) 

Using g ab = e~ x a ab and d t X = — v, we have 

^ ab = e- A (w a6 + ^a a6 ). (5.3) 

Combining these equations together with the Fuchsian expansion of a results 
in the equation: 

O t uf - e a u ab = 2e^~ £k)t a ac u b k:C , (5.4) 
which is of Fuchsian type if Sk > s a > 0. 

5.1.2 Equation for 

The Fuchsian expansion of k together with N = e x and d t \ = —v imply by 
straightforward computation that 

d t k b a ee e~ x { l -v 2 5 b a + e-^(v - e k )u b kA + e^ftu^}, (5.5) 
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and 

Nrk b a ee e-\ l -v 2 5 b a + ve'^u^) + e-^ k /-v5 b a + e-*«< a )}. (5.6) 

We see that e~ A -u 2 disappears from the difference d t k b a — Nrk b a , which moti- 
vates the choice of the Fuchsian expansion. 

To write the evolution equation 2.18 for k we now compute 

V b d a N = e- x a bc V c d a e x = a bc [d c Xd a X + d a d c \ - T d ac (g)d d \]. (5.7) 

On the other hand, since £ is 2 dimensional and g is conformal to a with a 
factor e A , we have that 

JVi£ = e x R b a ee ie^ifa) = i^{i2((7) - A CT A}. (5.8) 

From these results, if we define 

f a (t, u, u x ) := -V b d a N + NR b a - Nd a $.d b <5> 
then we calculate 

f b ee o- 6c [<9 c A<9 a A + d a d c \ - r^ c (^)9 d A] + l -5 b a [R{a) - A CT A] - a 6c $ a .$ c (5.9) 

We see that f b is at most a second order polynomial in t, is analytic in x 
when 5, w, A, 5" are analytic; is linear in du; and is analytic, bounded and 
Lipshitzian in u for u bounded and for large 4 t, except eventually for the last 
term which reads 

a bc <f> a .<f> c = 2a fec (0 a c + l -e-^u a u c ). (5.10) 

The expansion 4.10 of (f> a shows that it does not cause problems for the 
boundedness of f b . However the expansion of uj a gives 

e'^Ua = e~ 2(p d a (Co + e 2 ^) - e - 2 ^[2(l - cos9)d a + Q- l sin 2 6d a {® - wt)] 

+e-^' t u ula . (5.11) 

4 This restriction on t comes from the covariant components of a which remain bounded 
as long as a ab remains positive definite. 
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It follows from 4.6 that 

e 2(0-0) = with 6( p = e -e*t 

Therefore, using (1 — cos9)/sin9 = tan(9/2) we have: 

- e - 2 ^ - a _ _ 

e" 2 ^ = e 2 ^— — 9 a (w+e^)-e" 2 ^[2tan-9 a 0+e- 1 s2n^ a (0-^)]+e- £ "'V 

(5.12) 

We see that e 2( ^u; a will increase like (sinO) 1 - that is, like e^* - as t tends 
to infinity, except if 

uj + e 2 ^ = constant. (5.13) 

Condition 5.13 is a generalization of the condition imposed on the fields in 
[IM], with other notations, to obtain AVTD behaviour, in the case that S 
is a torus. Following the terminology of [IM] we call equation 5.13 the "half 
polarization" condition. Its geometric meaning is that the set of geodesies 
in the Poincare plane representing the VTD solution all tend to 
the same point of the axis Y = as t tends to infinity. 

After inserting the Fuchsian expansions and multiplying by e x+£k:t we find 
that the equation 2.18 takes the form 

&<« - e k u{ a - l -v5 b a ul c = e-^ul c u\, a + e x+ ^f b a (t, u, u x ). (5.14) 

Since A = A — vt and v — w this system can take a Fuchsian form only if the 
functions Co and satisfy the half polarization condition 5.13. To obtain the 
system in obviously Fuchsian form in that case, we split 5.14 into its trace 
and its traceless parts. For the trace part we have 

«, a - e k ul a - vul A = e-^X- + eX+£kt fa(t, u, u x ). (5.15) 

This equation takes Fuchsian form if and only if 5.13 is satisfied and v > e k . 
The same is verified for the traceless part T u b ka , which satisfies an equation 
with left hand side 

d t T u\ a -e k T u\ a . (5.16) 

5.1.3 Equation for u a >. 

Using the expansion of k and the relation d t X = —v, we find that 

d t a ab = 2e 2X k ab + a ab d t \ = 2e~ £k (5.17) 
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The equation for af' gives therefore the following equation for u c i : 

d t u«\ c - £ CT < C = 2e^'- £ '=) t V c K<c-] (5-18) 
which is of Fuchsian type so long as e a > < e k . 

5.2 Wave map equations. 

5.2.1 Equations for auxiliary variables. 

The equations resulting from the introduction of the new variables (f>t,u)t are 

^0-^ = 0, dt(v-u t = 0. (5.19) 
The first equation is of Fuchsian type for if £$ ( > £$ , since it reads 

d tU(t> - e^u* = e^+^u^. (5.20) 

The second equation reads 

[dtu u + (20 t - e u )uj[ - e<- e *- e **» Uut = 0. (5.21) 

We replace <p t by its value given in 4.8, which we write as follows 

<t> t = ~w + ^w(l - cosO) + e- £ **% t . (5.22) 

Since 1 — cosO falls off to zero as e~ 2wt , the equation 5.21 is of Fuchsian type 
for u u if w > and £$ ( > £$. 

In the equations 2.20 to be satisfied by </> a and u a , the derivatives of the 
VTD terms disappear, due to the commutation of partial derivatives. The 
equation for (f> a reads 

dtu^ - e^u^ = e- (e **- e *' ) *(9 o u 0t - td a w), (5.23) 

while the equation for u a becomes, using the expressions for u t and u a 

d tUu)a + (2<j) t - £u ,)u Ua = e-^- £ ^\d aUu)t + 2<j) aUiVt ). (5.24) 

These equations are of Fuchsian type so long as w > and £$ t > £$'. 
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5.2.2 Equation for u^ t . 

The first equation, 2.5, for the wave map reads 

g^{V a d P (t>+ l -e-^d a udpu) = 

-e' 2 \,9 t 0i + ^ e - 4 W)+e-V a6 (V a fe + ie- 4 ^ a ^)+^r° /3 t = (5.25) 

Using the Fuchsian expansions for cf> t and uo t together with 9' = —w sin 9 
and the value given in section 5.1.2 for e 2 ^ _< ^we find that: 

d t (f>t + ^e' H uj t uJ t = 

e~ e ^\d t u^ - e^Ufr) - ^w 2 sin 2 9 + ^(e~ 2 ^wsin9 + e'^u^f 
On the other hand, using the expansions for a ab , (fi a and u a we find that: 

e-V a6 (V a fe + ie- 4 ^) = 

e A {a ab +5a ab ){V b dJ+e- £ ^V b u^ a + 1 (C^+e"^'^ ) (e-^^+e" 6 "'*^ ) } 
We recall that 

va^ = v 6 [<9 a + — e- 1 9 a (e - at)], 

while under the half polarization assumption 

oj + e 2 ^ = constant, (5.26) 
the product e~ 2(i> d a Qj is given by 

e - 2 *d a (b = -e- 2S,t> Q- l sin9d a {e - w). 
Finally we calculate 

g^V% ee ^e- 2X = -e- 2 ^ul a . (5.27) 
14 



Inserting these computations into the first wave map equation produces 
an equation of the form 

dtUfa - e^Ufa = C^fa (x, t, u, du) (5.28) 

which is of the Fuchsian type 5.1 (with /x > 0) so long as v > £$ t , and 
Ek > £$ t - 

Analogous computations show that the equation for u Ut is Fuchsian pre- 
suming these same inequalities hold. 

5.3 Results for evolution. 

As a consequence of the calculations above we have proven the following 
theorem. 

Theorem 5.1 There exist a collection of positive numbers {e a , e a >, Ek,£<i>, £$ t , 
such that, given analytic asymptotic data on E, A = {v = w, A, a, 0,0, u}, 
the Einstein - wave map evolution system written in first order form for the 
unknown U, which defines g, k, $ and auxiliary variables by the Fuchsian 
expansions of section 4, is a Fuchsian system for U if and only if and Co 
satisfy the half polarisation condition 5.13 and v > 0. It admits then one and 
only one analytic solution tending to zero at infinity. 

To show that this result implies that we have a family of solutions of the 
Einstein - wave map evolution system which decays to solutions of the VTD 
equations, we need to verify that for a large enough t we have $t = <9i$, 
$ a = <9 a $ and the like. To show that $ a = (9 a $ we use the equations 2.20 
together with commutation of partial derivatives to show that: 

d t (<f>a - d a 4>) = d a <p t - d a d t cf> = 0; (5.29) 

hence <p a — d a (p is independent of t. As t tends to oo it tends to zero because 

0a - <9 a = e" £ ^'% a - e~ e **(0 o «0 - e^). (5.30) 

It must therefore always be zero. Analogous arguments can be used to show 
that u a = d a cu and of = V c a ah . 
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6 Constraints. 

The solution of the evolution system satisfies the full Einstein equations so 
long as it satisfies also the Einstein constraints, that is 

C : = £° = ~{R(g) - k.k + r 2 - e' 2X d t $ .d t ®} = 

C a := e A £° = -{V b k b - 8 a r + e'^.d^} = 0. 

As usual we will rely on the Bianchi identities, here to construct a Fuch- 
sian system satisfied by the constraints. Together with the wave equation 
satisfied by $, the Bianchi identities imply that 

(3) V a S« = 0. (6.1) 

Modulo the evolution equations ^R b a — p b a = that we have solved, with 
p b a ee $ a .$ 6 , it holds that 

®R-p = R°-p° ; (6.2) 

hence 

£° ee R° - pi - l -5^R - P )= l -5T ] R ~ P) (6-3) 

and 

K = -^R-p) = -5 b X (6.4) 

We use these equations and the identities 

S° ee e- x C a , H a ee -iV 2 S a0 ee -g ab N 2 Y? h = -e x g ab C b (6.5) 

together with the expressions for the Christoffel symbols of the metric ^g. 
We find that the equations 6.1 can be written in the form 

d t C - 2e x rC = g ab V a (e x C b ) + g ab e x d a \C a (6.6) 

and (after some simplifications and multiplying by e A ) 

d t C a - e x rC a = e x V a C + 2e x d a \C . (6.7) 

Equivalently, we have 

d t (e 2X C ) - 2ve 2X C - 2e x re 2X C = e x a ab V a (e x C b ) + a ab e x d a \e x C a (6.8) 
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and 

d t (e x C a ) - ve x C a - e 2X rC a = V a (e 2A C ). (6.9) 

We see that e 2A C and e x C a satisfy a linear homogeneous system, which 
admits zero as a solution. This solution is the unique one tending to zero at 
infinity, so long as the system is Fuchsian. 

Lemma 6.1 The system 6.8, 6.9 is Fuchsian, for a solution of the evolution 
system, if the VTD solution satisfies C = (i.e. v 2 = w 2 ). 

Proof. Since the coefficients of the equations 6.8, 6.9 are constructed 
from solutions of the evolution system we may use the expansions and esti- 
mates derived in previous sections. In particular we calculate 

e X T = v + e-^u a k!a . (6.10) 

Equation 6.8 can therefore be written as the following equation of Fuchsian 
type: 

d t (e 2X C ) - Ade 2X C e 2X C = e'^ut a e 2X C + e x a ab V a {e x C b ) + a ab e x d a \e x C a . 

(6.11) 

Equation 6.9 is not a priori in Fuchsian form for the pair (e x C a , e 2A C*o) in 
spite of the identity 6.10. However if we use the identity 

e 2X C = --{e 2X R(g) - e 2X k.k + e 2X r 2 - d t $.d t <Z>} (6.12) 
2 

and the property C = together with the expression for R(g) given in 5.9 
we can show that there exists a number fj, > and a bounded function F(x, t) 
such that we have 

|e 2A C | < e-**F(x,t) and |<9 a (e 2A C )| < e~ ltt F{x, t). (6.13) 

It follows that 6.9 takes Fuchsian form. ■ 

Theorem 6.2 A solution of the evolution system satisfies the full Einstein 
wave map equations if and only if the half polarized asymptotic data satisfies 
the condition w — v, and also 

= 1 and ve~ x+2 ^ = constant. (6.14) 
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Proof. To complete the proof that Cq = C a = it suffices to show that 
e 2A C*o and e x C a tend to zero at infinity. We have already checked that this 
is true for e 2A Co, as long as w = v; i.e. Co = 0. 

We now study the asymptotic behaviour of e x C a . If we denote by 5u the 
difference between a field u and its VTD value, we calculate (recall that 
A = A, a = a) 

e\C a - C a ) = e A {(V 6 - V b )k h a + V b 5k b a - djr} + <f(*A) (6.15) 

with 

8($ t .$ a ) ee 2<M0a + 20 a 50i + ^e- 2 ^5(e- 2 ^ a ) + e-^ooJie-^uJt). (6.16) 

We see that, in the half polarized case, the Fuchsian expansions imply that 
e x {C a — C a ) tends to zero as t tends to infinity. Using the expressions for k a 
and A, we see that e x C a reads: 

e x C a ee ±e x d a (e- x v) - $ t A = ±(d a v - vd a \ + vd a vt) - <Ma 

Using the expressions of A, $ t , $ a and the half polarization condition, we find 
after some computation that 

<Ma = -w{cos6dJ + ier 1 ^ - wt)}. (6.17) 

Thus we find that the terms containing t disappear from e x C a if v = w and 
6 = 1. It follows that e x C a tends to zero as t tends to infinity (recall that 
cosO tends then to 1) if and only if 

^[d a v -vd a X] -vd a <j) = 0, 
a condition equivalent to the hypothesis 6.14 of the theorem. ■ 

Remark 6.3 In the half polarized case the VTD solution only satisfies asymp- 
totically the VTD momentum constraint, and only after being multiplied e x . 
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